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^ ' Abstract 
O 

Q\ . Exact diagonalisations of the extended Hubbard model are performed. In the insulating 
regime it is shown that the nearest neighbour copper-oxygen repulsion, V, leads to Frenkel 
excitons in the charge transfer gap at values of V of the order of copper-oxygen hybridisation, 

O " 

t. In the metallic regime it is shown that the static charge-transfer and density- density 
correlation functions diverge as a function of V, indicating a charge-tansfer instability and 



phase separation. This is accompanied by a softening of the q — * mode of the dynamic 
correlation functions which is associated with the excitonic excitations responsible for the 
superconducting correlations observed in the proximity of the phase separation boundary of 
ref . (ITS . 
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The charge-transfer fluctuations Cu 2+ — > Cw + and 2 ~ — > 0~ coupled to the nearest 
neighbour copper-oxygen repulsion, V, leads to a variety of interesting phenomena. For 
example, in the insulating-stochiometric limit it leads to the formation of charge-transfer 
excitons [|TJ, ||, for which there is some experimental evidence from Raman || and optical 
@] experiments. In the metallic regime the softening of the charge-transfer excitations leads 
to a charge-transfer instability and phase separation, which suggests the possibility of an 
excitonic pairing mechanism ||, ||. This instability is also related to the pinning of the 
Fermi energy on hole doping, and to the creation of new states in the charge-transfer gap 
0. Such behaviour is qualitatively similar to the observations from photoemission || and 
O Is spectroscopy ||. 

The three band extended Hubbard model, which best describes this behaviour, has been 
studed by a number of authors. Littlewood et al. [JTtJ, following the initial suggestion of 



Varma et al. B, studied the charge-transfer excitations in the weak coupling limit, via the 



RPA. This work was followed up by Bang et al. |11J in the weak coupling limit, and by 



Raimondi et al. 0] in the strong coupling limit. These analyses all indicated superconduct- 
ing correlations, mediated by excitonic excitations, close to the charge-transfer instability. 
Barford and Long [DJ performed a unitary transformation to second order in the copper- 



oxygen hybridisation and found superconductivity and phase separation. Finally, Sano and 



Ono JT4|, Vermeulen et al. \TE\ and Stechel et al. [[Uj used the mapping onto a Luttinger 
liquid to show that the one dimensional model exhibits superconducting correlations and 
anomalous flux quantisation in the proximity of the phase separation boundary. 

It is the purpose of this letter to study in detail the dynamics of the charge-transfer 
excitations as a function of doping and V by solving exactly finite size clusters in one di- 



mension. We will study the charge-transfer excitons in the insulating regime. By performing 
finite size calculations we will compare the exact results to the predictions of the variational 
theory of Vermeulen and Barford 0. Next, by studying the dynamic charge-transfer and 
density-density correlation functions in the metallic regime we will relate the softening of the 
q — > mode and the divergence of the static correlation functions to the superconducting 
correlations and phase separation instabilities of the copper-oxide chain phase diagram of 



ref. 15 . 



Our model for the copper-oxide chain consists of two atoms per unit cell. Neglecting 
the oxygen-oxygen hybridisation and considering the Coulomb interaction up to first-nearest 
neighbours, the copper-oxide chain is described by the two band model Hamiltonian 



H = -t ( d l*Pjv + h.c.) + — XXpJ CT Pia - d\ a d ia ) 

<ij>a ija 

+ u dJ2 d h d ti4i d ii + U pJ2PnPnp]iPji + v d \A*p]a>Pjv'i (!) 

i j <ij>aa' 



where i and j are copper and oxygen sites respectively, < ij > represents nearest neighbours 
and the operator d\ a (Pj a ) creates a Cu (O) hole with spin a. A is the charge-transfer energy, 
Ud (Up) is the copper (oxygen) Coulomb repulsion and V and t are the copper-oxide Coulomb 
repulsion and hybridisation, respectively. We will denote N p as the number of particles and 
N s as the number of sites. Our choice of model parameters will be those relevant to the high 
temperature superconductors. Thus we take Ud = 8t, A = 2 — At and U p = — 4t where 



t ~ 1.5eV |fL7| . However, as we are interested in the role played by the nearest neighbour 
interaction, V, this will be left free parameter. 

Firstly we examine the model in the insulating phase at a hole density of n = 0.5. At 



this density the system is a charge transfer insulator for all A ^ 0, and the ground state 
is a spin density wave. The majority of the charge resides on the copper sites, forming a 
Neel state. The exciton is identified as the lowest frequency peak, E exc , of the dynamical 
current-current correlation function, J(w), when it appears inside the charge transfer gap, 
Egap- J(w) is defined by 

J(w) = --Im(G R (w)), 

71 

where 



G R (w) = F.T. < ip | j f (t)j(0) | Vq > 

H - (E + w) + ir] 



<^o\j 1 (-Tr 77^ )j\i>o>v->o- (2) 



G R (w) is the retarded Green function and j is the current operator defined as j = —iJ2ia 
(cj^q+io- — c\ +la ci a ). {| ip >, E } are the ground state eigensolutions. In practice we choose 
a finite value of r] to broaden the peaks. The Green function is calculated via the continued 
fraction technique The energy of the gap is defined by E gap = E(N p + 1) + E(N P — 



1) — 2E(N p ). Figure 1 shows a typical profile for the excitonic spectrum using a 12 site 
chain, and the values of Ud = 8t, U p — 0, A = 2t and V = 2t. The values of V cr u, the 
point at which E gap = E exc , are shown in tables (1) and (2) for 8, 12 and 16 site chains for 
A = 2t and 4t, respectively. The values for V cr u at iV = oo are found by finite size scaling, 
whereby Eg^ p (V) and E^ X S C (V) are plotted against and the straights lines are extrapolated 
t° = 0. Notice that V crit , which is of the order t, is a decreasing function of A and U p . 



The exciton may be pictured as a hole that has been excited from the lower Hubbard 
band (or conduction band), which is predominately copper in character, to the valence band, 
which is predominately oxygen in character. The real space picture is therefore of a hole on 
a copper site hopping onto a neighbouring oxygen site. This effectively leaves an 'electron' 
on the copper site and a hole on the oxygen site, and due to the Coulomb repulsion between 
neighbouring sites this leads to an effective attraction between the electron and the hole. 
This attraction reduces the energy to a point below the charge transfer gap and causes the 
formation of a tightly bound 'electron-hole' pair, or a Frenkel exciton. The energy in forming 
the exciton is a balance between the Coulomb attraction of the 'electron-hole' pair and the 
kinetic energy loss. The former is driven by V, whereas the latter is determined by t. Hence, 
when V is of the order of t we expect excitons to exist within the charge transfer gap. 

As a more quantitative study, we have extended the variational estimate by Vermeulen 
and Barford of V crit for the two dimensional copper-oxide plane. This was performed 
via a canonical transformation up to 0(t 2 ) in the hopping t of eqn(l), assuming strong 
coupling (Ud — > oo). Figure 2 shows V cr t t against A for the one dimensional infinite chain 
and several values of U p (where V = j, A = j and t = Notice that for large A (where 
the calculation becomes asymtopically rigorous) V cr u decreases as a function of U p , and is a 
constant function of A, i.e. V cr u decreases with A. This is in agreement with the numerical 
results of tables (1) and (2), and arises from the fact that increasing the oxygen repulsion 
increases the single-particle gap, but has little effect on the exciton energy. In contrast, the 
variational estimate of refpfl indicates that in two dimensions U p has no effect on V cr u for 
large A. We also confirm the validity of the variational estimate by comparing V cr u with the 
finite size numerical results in table (2), which indicates good agreement. 



We now turn to the role that the charge transfer excitations may play in the doped phase 



of the cuprate superconductors. From previous work [15| this model is well known to undergo 
phase separation within a hole density range of 0.5 < n < 1.0 and large enough V. There is 
also strong evidence that the system becomes superconducting in the proximity of the phase 
separation boundary. We intend to investigate the role played by correlated particle-hole 
excitations near the Fermi surface in these two physical effects. 

It is difficult to identify excitonic peaks from the current-current correlation function in 
the metallic regime. Instead, it is more convenient to study the dynamic density- density 
and charge-transfer correlation functions (i.e. susceptibilities) at low momenta. The dy- 
namic charge-transfer correlation function, in particular, focuses on local intra-cell charge 
fluctuations of copper and its surrounding oxygens. Such a local fluctuation in charge can 
be considered excitonic in character, and at low energy and momenta these will correspond 
to excitations near the Fermi surface. The dynamic charge-transfer correlation function is 
defined as 

X 5 g (w) = F.T. < 5l(t)5 q (0) >, (3) 

where 

= (4) 

V iV s I 

The dynamic density-density correlation function, Xq( w )i is defined by replacing 5 q by n q = 
J2k c l+q c k- I n terms of our lattice Hamiltonian n q is defined as -^^J2ii n i~ < n >)e~ tqXl . 
This function examines charge fluctuations across the whole lattice. From a thermodynamical 
argument it can be shown that the compressibility, k, is related to the static density- density 



correlation function in the long wavelength limit, , by 



C" o is the integrated dynamic correlation function (= J (a;)du;), p is the density, /? is 
the inverse temperature and V is the volume. So a divergence of the static density- density 
correlation function in the limit q — > is directly related to a divergence of k, and hence 
to phase separation. Moreover, the divergence of the static correlation function should be 
accompanied by a soft mode. Figures 3(a) and (b) show the position of the lowest peak in 
the dynamic charge-transfer correlation function and the integrated weight for q = jf- (i.e. 
the finite size q — > limit) on a 12 site chain at a filling of 8/12 and 10/12, respectively. As 
the nearest neighbour repulsion is increased the static charge-transfer correlation function 
rapidly increases and diverges at about V ~ 2.0t for 8 holes and V ~ l.bt for 10 holes. (The 
divergence is taken to be at the point of inflection, as a finite size system cannot have a 
truly diverging susceptibility.) This is accompanied by the softening of the diverging q — ► 0, 
uo — > mode. (All other q modes stiffen as V is increased.) Likewise, the static density- 
density correlation function also diverges. The frequency of the density-density soft mode 
is identical to the charge-transfer soft mode, as charge-transfer and density fluctuations are 
coupled. This arises from the strong Coulomb interactions on copper which means that 
local charge fluctuations will tend not to conserve local charge and so these excitations are 
coupled to non-local charge fluctuations. The divergence of the static correlation functions is 
consistent with the divergence of the discrete compressibility at V ~ 2.0t and 1M for 8 and 
10 holes, respectively, as shown in the phase diagram of ref . ITS] . Moreover, it is reasonable 



to associate the softening of the g — ^ charge-transfer excitations with the indications of 
superconducting correlations from both the Luttinger liquid charge exponent, K p , exceeding 



unity and to the onset of anomalous flux quantisation fl"5 |. 

The effect of U p is to suppress strongly both the static density-density and charge-transfer 
correlation functions. At n = A and V — 3t C™ drops from 0.22 at U p — to 0.09 at 
U p = 2t. 

In conclusion, we have shown that the nearest neighbour repulsion leads naturally to 
Frenkel excitons in the charge-transfer gap at values of V cr u of the order of t in the insu- 
lating regime. Further, we have found that in one dimension increasing U p has the effect of 
decreasing the value of V cr u, and we have given a qualitative argument to explain this effect. 
These numerical results are in good agreement with the analytical results of 0. 

In the metallic regime we have shown that the static charge-transfer and density- density 
correlation functions diverge as a function of V, indicating a charge-tansfer instability and 
phase separation. This is accompanied by a softening of the q — > mode in the dynamic 
correlation functions which we associate with the excitonic excitations responsible for the 
superconducting correlations observed in the proximity of the phase separation boundary of 
ref . [115 . 
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Table 1. 



A = 2* 
U d = 8t 


Number of lattice site, N s 


8 


12 


16 




oo 


u p = o 


2.06* 


1.7t 


1.54* 




1.29* 


U p = At 


1.57t 


1.22* 


1.09* 




0.88* 



Table 2. 



A 


= At 


Number of lattice site, N s 


u d 


= 8* 


8 


12 


16 




oo 


Analytic 


u p 


= 


2.02* 


1.54* 


1.38* 




0.94* 


1.22* 


Up 


= At 


1.82* 


1.31* 


1.16* 




0.77* 


1.05* 



Figure Captions 



Figure 1. The current-current response function at Ud = 8t, U p = Ot, A = It and V — It 
for the stoichiometric 12 site chain. E g is the charge transfer gap energy and E e identifies 
the exciton energy. 

Figure 2. V cr i t (= V cr i t /t) versus A(= A/i), where t = t 2 /A, for the one dimensional chain, 
using the theory of ref. ||, in the Ud = oo limit. The exciton is a bound state above the 
lines. 

Figure 3. (a) The q = ? static charge-transfer (diamonds) and density-density (solid 
squares) correlation functions, and the position of the lowest frequency mode of the dynamic 
correlation functions (open squares) as a function of V/t for 8 holes in a 12 site chain. The 
frequency is in units of t. Notice the different normalisations for the correlation functions. 
Ud = 9t, Up = and A = 2t. (b) The same as (a) with 10 holes. 

Table 1. Values of V cr u for various system sizes. 

Table 2. Values of V cr u for various system sizes. The analytic result is in the Ud = oo limit. 



